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ABSTRACT 

Monte Carlo simulation is one of the best tools to study the complex spectra 
of Compton-thick AGNs and to figure out the relation between their nuclear 
structures and X-ray spectra. We have simulated X-ray spectra of Compton- 
thick AGNs obscured by an accretion torus whose structure is characterized by a 
half-opening angle, an inclination angle of the torus relative to the observer, and 
a column density along the equatorial plane. We divided the simulated spectra 
into three components: one direct component, an absorbed reflection component 
and an unabsorbed reflection component. We then deduced the dependencies 
of these components on the parameters describing the structure of the torus. 
Our simulation results were applied to fit the wide-band spectrum of the Seyfert 
2 galaxy Mrk 3 obtained by Suzaku. The spectral analysis indicates that we 
observe the nucleus along a line of sight intercepting the torus near its edge, and 
the column density along the equatorial plane was estimated to be ~10^^ cm~^. 
Using this model, we can estimate the luminosities of both the direct emission 
and the emission irradiating the surrounding matter. This is useful to find the 
time variability and time lag between the direct and reflected light. 

Subject headings: galaxies: active — galaxies: Seyfert — radiative transfer — 
X-rays: individual (Mrk3) 



1. Introduction 

Active galactic nuclei (AGN) emit huge amounts of energy as a result of accretion of 
matter onto supermassive black holes (hereafter SMBHs). Strong hard X-ray emission from 
these nuclei can be represented by power-law emission with a canonical photon index of ~1.9. 
The strong emission irradiates material around the SMBH, leading to reprocessing of the 
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emission. Observation of the reprocessed emission is conducive to revealing the environment 
of the black hole, and hence help to understand the problems of the fuel supply and evolution 
of SMBHs. 

Compton-thick AGNs are suitable for the study of the reprocessed emission in AGNs, 
since their reprocessed emission dominates over the direct emission below 10 keV due to 
large column density A^h > 10^^ cm~ ^. In addition, Com pton-thick AGNs are thought to be 
abundant in the local universe (e.g.. iRisaliti et al.lll999l ). Thus, these objects are important 
for understanding some key problem s in AGN researc h, for example synthesis modeling of 
the Cosmic X-ray Background (e.g., lUeda et al.ll2003l ) and the growth of AGNs. However, 



their spectra are so complex that the detailed nature of Compton-thick AGNs has thus 
far been unclear. There is another problem for spectral analyses of Compton-thick AGNs. 
The baseline model based on previous observations, which consists of both a direct and a 
reflection compon ent, has been used to reproduce the complex X -ray spectra of Compton- 
thick AGNs (e.g.. iMagdziarz fc Zdziarski Ill995l : ICappi et al.lll999f ). This model worked well 
for reproducing their spectra, and has contributed to our understanding of AGN. However, 
it is difficult to obtain information about the structure of the surrounding material from 
spectral fitting with this baseline model because the reflection model was developed for 
an accretion disc geometry. Thus, the reflection model does not match exactly with the 
actual reflection from the material around the black hole. Therefore we need a model that 
represents actual X-ray emission from the surrounding material using detailed Monte Carlo 
simulations. 

Previous work that has investig ated AGN X-ray spectra by means of Monte Carlo sim- 
ulations so far includes the following. iGeorge fc Fabian! (jl99ll ) assumed a semi- inflnite, plane- 
parall el configuration, and they calculated refiection from a plane-parallel slab. IWilman fc Fabian 
(Il999l ) computed transmission thr ough a homogeneou s sp here of cold mate r ial wi th a radius 
determined by its column density. lAwaki et al.l (Il99ll ) and lGhisellini et al.l (119941 ) computed 
X-ray spectra emitted from AGNs and assumed a torus structure which has an open refiect- 
ing area characterized by its half-opening angle. However, these results were not applied to 
make available convenient spectral-fitting models to reproduce the X-ray spectrum of, for 
example, a Seyfert 2 galaxy. 

We assume surrounding material with a 3D torus configuration and then simulate a 
spectrum from an AGN, considering the effect of Compton down-scattering and absorption. 
Performing Monte Carlo simulations, we investigate the relation between continuum emis- 
sions and the structure of the torus. Furthern iore, the strengjth of iron line is important to 



reveal the existence of a hidden nucleus (e.g., iMaiolino et al.lll998l ). Thus, we also investi- 



gate a dependency of the iron line on the structure of the torus. We apply our simulation 
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results to the spectral fitting of the X-ray spectrum of Mrk 3 observed with Suzaku, thereby 
constraining the structure of the torus in this Seyfert 2 galaxy. 



Model Definition and Calculation 



2.1. Basic Assumptions 

In our Monte Carlo simulation, we adopted a standard spherical coordinate system 
with the radial distance (r) from the origin, the zenith angle {6) from the Z-axis, and the 
azimuthal angle (0) from the X-axis. A primary radiation source was placed at the origin, 
and illuminates the surrounding material. The material was assumed to be neutral and cold 
(T < 10" K). 

We took into account photoelectric absorption, iron fluorescence, and Compton scatter- 
ing in our simulation. The photoelectric absorption cross-section, crabs, w as calculated by us- 
i ng th e NIST XCOM database and the cosmic elemental abundances of Anders fc Grevesse 



( 119891 ). We also calculated the photoelectric abso rption cross-section of Fe for deciding iron 
absorption events. Note that the cross-section by iBalucinska-Church fc McCammon I (119921 ) 
is not valid for energies ab ove 10 keV, an d that the NIST cross-section in the 1-100 keV band 
is nearly equal to that by IVerner et al. I (119961 ) which is identified as vern in XSPEC (see 
Figure 1). The Compton scattering cross-section, cTes, was calculated with the Klein-Nishina 
formula. The number density of electrons, rigs, for Compton scattering was related to the 
effective hydrogen number density, nn, by n^s = 1.2 nn. 



We used the iron K-shell fluorescence yield of 0.34 (IBambynek et al.lll972h and a rati o 
17:150 between the iron and Kq, fluorescence line transition probabilities (jKikoirull976l ). 
Although the iron Kq, fluorescence line consists of two components, K^i and K^y, at 6.404 



and 6 .391 keV, respectively (for neutral iron) with a branching ratio of 2:1 (IBambynek et al. 
19721 ). we made no distinction between Kd and Ka2 photons, adopting a common value of 
6.40 keV. The iron Kp line is 7.06 keV for neutral iron. 



2.2. Monte Carlo Simulation 



A cross section view of the adopted structure of the torus is illustrated in Figure 2. The 
center of the torus is placed at the origin of the coordinate system, and the equatorial plane 



ihttp://physics.mst.gov/PhysRcfData/Xcom/Text/XCOM.html 
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of the torus structure lies in the X-Y plane. The structure of the torus is described by the 
following structure parameters: the half-opening angle 6'oa, the column density A^h along the 
equatorial plane, and inner (rjn) and outer (rout) radii of the torus. We assumed rin/rout = 
0.01 in our simulation. The angle di in Figure 2 presents the inclination angle of the torus 
relative to the observer. 

A ray-trace method was adopted in our Monte Carlo simulation. The primary X-ray 
source was assumed to emit photons with the energy spectrum, I{E) oc i?~°'^exp(— i?/i?c), 
where Ef., is a cut-o ff energy fixed at 360 keV. The spectrum is typical for Type 1 AGNs (e.g.. 



Madau et al.lll994l ). The primary X-ray source was assumed to be isotropic. Each photon 
from the primary source had both an initial energy and an initial direction of propagation. 
In the case that a photon was injected into the torus, an interaction point of the photon was 
calculated by using a random number (see below). If the Compton scattering occurred at 
that point, the energy and the direction of photon was changed. The photon was tracked 
until it escaped the torus structure, or until it was absorbed in the torus. Note that when 
the photon with an energy above iron K-edge was absorbed by iron, a K-shell fluorescence 
line was isotropically emitted with a probability of the K-shell fluorescence yield in our 
simulation, where the iron absorption event was decided by using the ratio between the 
photoelectric absorption of iron and (Tabs- Both energies and directions of propagation of all 
escaping photons were recorded in a photon list. For a given observed torus inclination angle 
of 9i (see Figure 2), we extracted photons whose zenith angles of propagation ranged within 
6'i± 1°, from the photon list, and then the extracted photons were accumulated into energy 
bins to form a spectrum. 

The photon transportation in the torus is a key technique in our Monte Carlo simulation. 
The distance / to the next interaction is determined by the probability p, which is described 
as follows: 



p = exp(-r) = / exp(-atot nn)dl, (1) 
J 

where r and Utot are an optical depth and the total cross section of the interaction, re- 
spectively. The (Jtot is comprised of the sum of cxabs and aes- By inverting the cumulative 
probability function, / is expressed as 

Z = ^— = ^xln(p). (2) 

nno-tot nno-tot 

The distance / is calculated from a uniform random number between and 1 referred to as 
p. 
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Another key in our simulation is Compton scattering. The scattering angle, ^scat (relative 
to its direction of propagation), is calculated by the differential cross-section for Compton 
scattering. We assumed that the differential cross-section was proportional to (l+cos^ ^'scat), 
as for the Thomson differential cross-section, since th is approximation is effic ient for analysis 
of the Suzaku data below a few hundreds keV (e.g., iGeorge fc Fabianlll99l[ ). We note that 
this assumption over-estimates back-scattering relative to forward scattering at high energies. 
The effect of this approximation is seen in the reflection components, and depends on the 
geometry and A^h of the torus. In order to estimate the effect, we performed a simulation in 
the case of ^oa=40°, 6'i=41°, and A^h=10^^ cm~^, and found that the change of the reflection 
components due to this assumption was ~10% at 100 keV. 

The azimuthal angle of the scattering, (pscat, was randomly selected in the region 0° < 
0scat < 360°. The energy of the scattered photon was changed to be Ei^/ {1+^^(1— cos9 scut 
)), where m and are electron mass and en ergy of the photon before Compton scattering, 
respectively (e.g.. iRybicki fc Lightmanlll979l ). 



3. Results 

3.1. Simulation Results 

Figure 3 is an example of a simulated AGN spectrum with Ah = 10^^ cm~^, 6'oa = 40°, 
and 6i = 45°. We generated a total number of 2.5 x 10* photons. The same number of photons 
were generated for each run throughout this paper. We separated the simulated spectrum 
into direct and reflection components, where the direct component has no interaction with 
the surrounding material, while the reflection component consists of X-ray photons reflected 
in the surrounding material. Since an reflection component, which was modeled by pe xrav 



in XSPEC, was required in the Suzaku spectrum of Mrk 3 by lAwaki et al.l (120081 ). we 
divided the simulated reflection component into two, which are referred to as the reflection 
components 1 and 2. The reflection component 2 consists of photons emitted from the inner 
wall of the torus without obscuration by the torus. The reflection component 1 consists of 
the rest of the reflected photons (see Figure 2). We have studied the parameter dependence 
of these three components based on the simulations with various Ah? ^oa? and d[. 
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3.2. Dependence of the Continuum Emission on the Structure Parameters 

3.2.1. Nn Dependence 

For studying the A^h dependence, we simulated spectra of the three components with A^h 
= 5xl02^,102'',2xl02'',3xl02'',5xl02^, and 10^^ cm^^ and the simulated components are 
shown in Figure 4. In these runs, we set ^oa = 40° and 9i — 50°. All the three components, 
especially the direct component, show a dependence on A^h- It is expected that the direct 
component is affected only by the column density (A^^*^) along our line of sight. We compared 
the simulated direct components with the cut-off power-law models, which were affected by 
both photoelectric absorption and Compton scattering (Figure 5) . It is found that the models 
were in good agreement with the simulated direct components. Please note that the column 
density N^^ is a function of 9i, ^oa and the ratio r-m/rout (=^)- The ratio of A^^^ to the A^h is 
described as 

j^is,j^ _ rjcos 0i - cos gpa) + sin(gi - gpa) , . 

In the case of r=0.01, 0oa = 40°, and 0i = 50°, for example, the value of N^^/Nji is deduced 
to be 0.97. 



3.2.2. Dependence on the Half-opening Angle, 0os, 

Three components depended on the half-opening angle, 9o^, but we discuss the 9oa- 
dependence of only the reffection components 1 and 2, since the ^^oa-dependence of the direct 
component, which is affected by A^'*, has been described by equation 3. 

The left panel in Figure 6 shows simulated spectra of the reflection component 1 for ^oa 
= 10° — 70° in steps of 10°. The values of A^h and 9\ were fixed at 10^^ cm~^ and 90° in these 
runs, respectively. The spectral shape of the reflection component 1 has a small dependence 
on 6*03, and the total photon count in the reflection component 1 decreases with increasing 
^oa- We plot the total count in Figure 7. Since some photons emitted from the central source 
escape through the opening area of the torus, the intensity of the reflection component 1 
is dependent upon ^oa- We compared the total counts of this component with the covering 
fraction of the torus, cos^oa- In the right panel in Figure 8, we plot the normalized total 
counts, which are divided by the total counts in the case of spherical distribution of the 
surrounding material. It is found that the ^^oa-dependence of the reflection component 1 can 
be roughly explained by the covering fraction of the surrounding material. For large A^h 
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of 10^^ cm ^, the total counts do not follow this ^oa-dependence as for low A^h due to the 
absorption of the reflection component by the torus itself. 

The right panel in Figure 6 shows the ^oa-dependence of the reflection component 2. In 
these runs, we set ^i=^oa+l° in order to obtain a high intensity of the reflection component 
2 under the condition that the nucleus is obscured by the torus. The total count of the 
reflection component 2 are shown in Figure 7. It is found that the reflection component 2 is 
to be zero at both Ooe, =0 and 90°, and has a maximum intensity at ^^oa ~ 30°. The reflection 
component 2 is emitted from the visible inner wall of the torus. Thus, this component 
vanishes at ^^oa =0 due to there being no inner wall. This component is also associated with 
the total number of photons injected into the torus from the central source, which is related 
to the solid angle of the torus. Therefore, the intensity of the reflection component 2 is 
associated with a combination of both apparent size of the visible inner wall of the torus 
and sohd angle of the torus, i.e. cos^oacos(2^oa — 7r/2). The left panel in Figure 8 shows the 
total counts of the reflection component 2 in A^h=10^^, 10^^, and 10^^ cm~^. For the case 
of A^H=10^^ cm~^, the ^oa- dependence is well reproduced by the combination described as 
cos^oacos(2^oa ~ T^/'^)-i sincc the reflection occurs mainly on the inner surface of the torus. 



3.2.3. Dependence on the Inclination Angle 9\ 

We show simulated spectra of the reflection components for ^oa=10° with to 
71° in steps of 20° in Figure 9, and show the total counts at ^oa=10° and 30° against 
from 1° to 89° in steps of 2° in Figure 10. In these runs, A^h was fixed at 10^^ cm~^. The 
reflection component 1 has a weak ^^j- dependence (see Figure 9), and has a similar intensity 
in > ^oa +10° (Figure 10). The ^^j- dependence is consistent with the fact that the reflection 
component 1 depends on the covering factor of the surrounding material as mentioned in the 
previous subsection. The unified model of Seyfert galaxies predicts that 6i is larger than 6oa 
in Seyfert 2 galaxies. Thus, for most Seyfert 2 galaxies, the reflection component 1 should 
show a small ^i-dependence, and this component should depend on mainly ^oa ^-nd A^h of the 
torus. In the case of ^oa < < ^oa+10°, the reflection component 1 should have a strong 
^i-dependence. 

The reflection component 2 shows a strong ^i-dependence (Figure 9) , because the com- 
ponent is emitted from the visible inner wall, whose area is decreasing with increasing 9i. If 
we observe the source edge-on, i.e.; 6\ =90°, the area of the wall is apparently zero for the 
observer. On the other hand, if we observe the source face-on, i.e.; ^i=0°, the projected area 
(i.e. the area projected onto the sky, perpendicular to the line of sight) is largest. The 9\- 
dependence of the reflection component 2 is roughly explained only by that of the projected 
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area. Note that larger portion of the projected area, or the visible inner region, is obscured 
by the near side of the torus with increasing 9i in 9i > Ooa- The spectral shape of the reflection 
component 2 is also affected by 6i as shown in Figure 9. With increasing 9[, the low-energy 
cut off becomes more pronounced. This is because only X-rays reflected from the visible 
inner region of the torus (see Figure 2) contribute to the reflection component 2, as dcflncd. 
Some photons emitted from the central source enter the far side of the torus and Compton 
scattered. A part of the scattered photons are absorbed inside the torus before they ever 
reach the visible inner region. Since the length of the path of such photons increases with 
9i, more photons are absorbed before they reach the visible region. As a result, we can see 
a ^i-dependence on the spectral shape. 



3.2.4- Note on dependence on rin/rout 

The ratio Tin/ rout (=^) was fixed at 0.01 in our simulation. We briefly discuss the spectral 
dependence on this ratio. We simulated the reflection components with different values of 
the ratio, 0.001, 0.01, 0.1, and 1. In these runs, A^h, ^oa, and 9i were set to be 10^^ cm~^, 40°, 
and 41°, respectively. The reflection component 1 shows a very weak dependence on this 
ratio, since these simulations were carried out for the same 9os., and this component depends 
mainly on the covering fraction of the torus. On the other hand, the reflection component 2 
below ~ 4 keV depends on the ratio r (Figure 11). Reflections occur more frequently than 
absorption at the inner region of the torus. Thus, in the case of the small r, the closest region 
to the source is obscured by the torus, since we consider the case of Seyfert 2 galaxies, 9i > 
9o3.. This effect is also seen in the 9i dependence of the reflection component 2 (see Figure 9). 
This study indicates that it is difficult to produce a strong unabsorbed refiection-component 
in our simple torus geometry with a small r (=0.001). 



3.3. X-ray luminosity absorbed by the torus 

In our simulation. X-rays emitted from the central source were absorbed by the dusty 
torus. It is important to compare the absorbed X-ray luminosity with the infrared lumi- 
nosity, since the infrared luminosity is a good indicator for the absorbing luminosity in the 
torus. We estimated the fraction of the absorbed luminosity with respect to the intrinsic 
source luminosity. The intrinsic luminosity of the central X-ray source was calculated from 
the assumed spectrum in section 2.2, and the absorbed X-ray luminosity was derived from 
the subtraction of an output luminosity from the intrinsic luminosity, where the output lu- 
minosity was deduced from a total energy of the escaping photons from the torus structure. 
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Figure 12 shows the absorbing fraction to the 1-100 keV intrinsic luminosity as a function 
of ^oa in A^H=10^^, 10^^, 10^^, and 10^^ cm~^. The dependence of the absorbing fraction is 
roughly explained by cos ^oa? which displays the covering fraction of the dusty torus. In the 
case of iVn = 1x10^^ cm~^ and 6'oa=45°, the absorbing fraction of 0.37 is obtained, and the 
absorbed luminosity is estimated to be 1.1x10^^ erg s~^ for the 1-100 keV intrinsic luminos- 
ity of 3x10^^ erg s~^, which corresponds to the 2-10 keV intrinsic luminosity of 1x10^^ erg 

The infrared luminosity has a good correlation with the 2-10 keV luminosity. By using 
the relation between these luminosities by Mulchaey et al. (1994), the infrared luminosity 
for an AGN with the 2-10 keV X-ray luminosity of 1x10"^^ erg s~^ is estimated to be 10"^^ 
erg s~^, which is about 10 times larger than the absorbed X-ray luminosity. The dust in 
the torus is heated by optical and UV photons as well as X-ray photons. Our study shows 
an estimation of the fraction of the absorbed X-ray luminosity to dust heating. We indicate 
that the absorbed luminosity depends on the geometry of the torus, 6'oa in Figure 12. Lutz et 
al. (2004) pointed out that the scatter in the relation between mid-infrared and absorption 
corrected hard X-ray luminosities was about one order of magnitude, and that the scatter 
was likely caused by the geometry of the absorbing dust. The scatter seen in the relation 
may be explained by the ^oa-dependence of the absorbed luminosity. 



3.4. Dependence of the Iron-line Emission on the Structure Parameters 



A prominent iron line with an equivalent width of > a few 100 eV is an important 
characteristic of Seyfert 2 galaxies. The dependence of the equivalent width of th e iron line on 



the s t ructure parameters of the torus have been studied in previous work (e.g.. lAwaki et al. 



199ll : iGhisellini et al.lll994l : iLevenson et al.ll2002l ). These studies were mainly performed for 
the equivalent width relative to the total continuum emission, comprised of the sum of the 
direct and reflection components. In our study, we have investigated the equivalent width 
relative to the reflection component (hereafter EW^ei) as well as the equivalent width to the 
total continuum emission (hereafter EWtot)- In Figures 13 and 14, we plot the equivalent 
width as a function of A^h and 6'oa respectively (note that for 6'oa=0 we use a spherical 
distribution instead of using our simple torus model with 6'oa=0). 

The left panel in Figure 13 shows the A^n-dependence of EWtot fixed at ^oa=30°. We 
found that the A^n-dependence is similar to those o f previous studi e s, alth ough our results 
are about 1.7 times larger than those obtained by IGhisellini et al.l (119941 ). The difference 
may be caused by the difference of iron abundance used in the simulations. The right panel 
in Figure 13 shows the 6'oa-dependence of EWtot- The EWtot in the Compton-thin region 
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( A^H < 10^^ cm~^) decreases with increasing ^oa- This ^oa-dependence is represented by 
a function of cos^^oa? similar to the reflection component 1, since the iron line intensity is 
proportional to the solid angle subtended by the surrounding matter at the source. 

The left panel in Figure 14 shows the A^n-dependence of EW^ei for 6*00=30°. The EWrei 
shows little A^h dependence in the A^h <10^^ cm~^, since the EWrei is mainly determined 
by the ratio between the absorption and scattering cross-sections. In the Compton-thick 
region, the EWj-ef shows large 6'oa- and 6'i-dependences. This is caused by the fact that the 
contribution of the reflection component 1 to the reflection continuum at the iron band is 
decreasing with increasing A^h (see Figure 4). 

We noted that the EW^ei in the Compton-thick regime is greater than 1000 eV in our 
simulations. Thus, observed lower values of EW^ei, less than 1000 eV, may indicate a low 
metal abundance of iron. 



4. Application to Observed Spectrum 

By means of Monte Carlo simulations, we showed that the three continuum components 
depend on the structure parameters of the torus, A^h, ^oa, and 9i. Our study suggests that we 
can estimate the structure of the torus by determining the three components in an observed 
spectrum. We therefore made a new model for spectral fitting based on our simulations 
with Ix 10^ photons in each run. In our new model, the direct component is reproduced by 
phabs * compcabs * {cutof fpl) in XSPEC, where compcabs is a new model that we made for 
representing Compton scattering, and the two reflection components are reproduced by using 
table models, which have parameters of photon index, A^h, ^oa, and 9i, since the reflection 
components are too complex to be represented by numerical expression. The table models 
cover the ranges of photon index of 1.5 — 2.5, A^h of 10^^ — 10^^ cm~^, 6'oa of 0° — 70°, and 
of — 90°. The details of the parameters of the table models are listed in Table 1. 



The Suzaku satellite can obtain a wide-band spectrum with good quality (IMitsuda et al. 



20071 ). We applied the new model to the wide-band spectrum of Mrk 3 observed by the 



Suzaku satellite on 2005 October 22-23 durin g the SWG p h ase. W e obtained XIS and HXD 



spectra in the same manner as described by lAwaki et al.l (120081 ). and then simultaneously 



fitted the XIS and HXD spectra above 1 keV with the new model defined as follows: 

/(ph s-i cm-2 keV^^) = PLl+e-'^^^Hig-^csAfHip2,2+re//ectionl(A'EL)+re//ection2(A/'H2)+^ivS, 

(4) 

where PLl and PL2 are power law components with a high energy cut-off. The high energy 
cutoffs of both PLl and PL2 were fixed at 360 keV, which is consistent with the cut-off energy 
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(Ecut) of the power law component Ecut > 200 keV obtained by ICappi et al.l (Il999l ). A^hi 
is the column density along our line of sight. We used the absor ption cross-section of vern 
for (Ja in the phabs model in XSPEC vl2.4. The abundances of lAnders fc Grevessd (119891 ) 
were used. The a^s is a Compton scattering cross-section, which is used in the compcabs 
model. The reflection components 1 and 2 in our simulation were reproduced by the two 
table models, reflectionl and reflection2, respectively. We set their photon indices equal 
to that of PL2, while we did not link their nomalizations to that of PL2, and their column 
densities along the equatorial plane {N112) were not linked to the column density (A^hi) of 
PL2 in our spectral fit. The emission lines (ELs) seen in the spectrum were represent by the 
sum of Gaussian components: 

ELs = gauss{Ei, o"i, A^i), (5) 

i 

where E^, cTj, and are the cent er energy, width, an d intensity of the i-th line. We fixed E^ 



and (jj at the values obtained by lAwaki et al.l (120081 ). 



We fitted the spectra with the new table models in the energy range from 1 to 70 keV. 
Since the table models were generated from simulations, the table models have statistical 
deviations. A typical standard deviation per 20 eV bin of the sum of reflection 1 and 2 in 
the 3-5 keV band is about 2%, which is about 1/5 of the statistical error of the observed 
data in this energy band. Due to the deviations of the table models, the "wiH have a 
fluctuation of about (bin number) x (|)^, which is estimated to be about 5. Thus, it is hard 
to find the best-fit parameters and their confidence regions with the x^-fitting procedure. 
We performed x^-test with our table models on the parameter grids of 60a, and 6i in order 
to examine whether we can obtain the structure parameters from the spectral fit with the 
table models. Table 2 lists the on the grids with 612 d.o.f. The minimum (Xmin) 
obtained for 9i < 6'oa, and the was more than Xmin + 30 for 9i > 6'oa+3°. Since the 6*1 
must be larger than the 60a, in our simple torus model (otherwise the Nm becomes zero), the 
study indicates that we observed the Mrk 3 nucleus near the edge of the torus. This is 
expected from the strong unabsorbed reflection component of Mrk 3 in the baseline model. 

On the other hand, we found that it is difficult to constrain the half-opening angle 60a 
from the spectral analysis, since the shapes of the reflection components 1 and 2 show a little 
6'oa-dependence (Figure 6). Although their intensities depend on ^oa as shown in Figure 7, the 
intensities of the reflection components also depend on the luminosity of the central source. 
Figure 15 shows the change of the normalization of the direct and reflection components for 
^i=^oa+l°- The normalization of the reflection component is roughly proportional 1/cos^oa- 
This relation indicates that the contribution of the reflection components to the observed 
spectrum is nearly constant. Furthermore, the spectral shape of the reflection component 
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1 in the 5-70 keV band is similar to that of the direct component at A^h ~10 cm~ (see 
Figure 3). As a resuh, it is hard to constrain ^oa? even if we hnk the intensities between the 
direct and reflection components. We note that the iron EW is not helpful to find the 6'oa of 
Mrk 3, since -EPKef has a little dependence on 6'oa in < 10^*^ cm~^. 

In order to constrain A^h2, we fitted the observed spectrum with our model on the grid 
of A^H2- Since we did not constrain the 6'oa from the spectral fit, we fixed ^oa at 50°, and 9i 
at ^oa+l°- The value of A^hi was a free parameter in this fit. We found that there is a 
minimum around ~10^^ cm~^, and that the is greater than x^jj^+30 for A^h2 < 6x10^^ 
and for A^h2 > 2x10^^ cm~^. We found that we constrain the structure parameter A^h of the 
torus by using our new model. 



The opening angle may be estimated from the opening angle of the NLR. ICapetti et al. 



(119951 ) found a NLR opening an gle in Mrk 3 of > 100°. The half opening angle of >50° 
is larger than the estimation by iRuiz et al.l (120011 ) d ue to the inclusion o f all the Z-shape 
emission components in the NLR in the estimation of ICapetti et al.l (119951 ). We here set 9oa 
= 50° and 6i = 51° in our spectral fit. Furthermore, the column density A^hi was linked with 
A^H2 by A^Hi=0.74 A^h2, by using the equation (3). We obtained a small reduced value 
of 1.18 (d.o.f.=613), which is comparable to that with the baseline model. Table 3 shows 
the best-fit parameters, and Figure 16 shows the best-fit spectrum. The photon index of 
the power law component and the column density of A^hi were estimated to be ~1.82 and 
~l.lxlO^'^ cm~^, respectively. The intrinsic luminosity of the power law component in the 
2-10 keV band was estimated to be 2.1x10'^^ erg s~^, which is about 1.3 times that of the 
estimate wit h the baseline mo del and is consistent with that obtained by using the model 
with plcabs (lAwaki et al.ll2008l], which describ es the X-ray transmission, correctly taking into 
account Compton scattering (lYaqoob 1119971 ). On the other hand, the intrinsic luminosity 



irradiating the accretion torus is estimated to be 5.1x10^^ erg s~^ from the normalization of 
the reflection component. The discrepancy of the intrinsic luminosities between the direct 
and reflection components may be arisen by a time lag of the reflection component, since 
the long time variability of Mrk 3 has been reported and t he distance of the accretion torus 
is estimated to be greater than 1 pc from the center (e.g.. 



Awaki et all l200d . 120081 ). 



5. Summary and Conclusion 

We simulated AGN spectra by using the ray-trace method, and made a new model 
for fitting spectra of Compton-thick AGNs. In our simulations we assumed an accretion 
torus surrounding a nucleus, which was characterized by 6'oa, 6'i, A^h, and the ratio of the 
inner and outer radii of the torus. We considered interactions of photoelectric absorption. 
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iron fluorescence, and Compton scattering in the simulation. The simulated spectra were 
separated into three components: one direct component and two reflection components. 

The direct component consists of X-ray photons which have no interaction with the 
surrounding material, and the component is only affected by column density along our line 
of sight. In fact, the direct component was well modeled with a cut-off power law emission 
affected by both photoelectric absorption and Compton scattering. On the other hand, the 
reflection components had not only an A^n-dependence but also both Oq.^- and ^^j-dependences. 
The reflection component 1 shows a ^^oa- dependence, which is explained by the covering factor 
of the torus. The reflection component 2 has dependence on both and 9\. The dependence 
of the reflection component 2 is roughly explained by the apparent size of the visible inner 
wall of the torus. 

We fltted the Suzaku Mrk 3 spectrum with the new model based on our simulations, 
and found that the spectrum could be represented by our model. The structure parameters 
of the torus of Mrk 3 were estimated with the new model: 9\ ~ ^oa + 1° and A^h ~ 10^^ cm~^, 
although it was hard to constrain ^oa from our spectral analysis. We estimated the intrinsic 
luminosity of the direct component and the intrinsic luminosity irradiating the surrounding 
matter. Assuming 9^^—bQ° and A^hi=0.76A^h2, the 2-10 keV luminosity of the direct compo- 
nent was estimated to be about 2/5 of that irradiating the surrounding matter. This may be 
explained by time variability of Mrk 3 and time lag between the direct and reflected lights. 

We demonstrated that we can bring out the structure of the torus from an observed X- 
ray spectrum with our new model. The wide-band X-ray spectra will be helpful to determine 
the structure of AGNs. 

We also thank Drs. T. Yaqoob and K. Murphy for useful discussions and careful reading. 
We also thank the anonymous referee for helpful comments and suggestions. This study is 
carried out in part by the Grant support for Scientiflc Research of Ehime university (H.A.) 
and the Grant-in- Aid for Scientiflc Research (17740124 Y.T.) of the Ministry of Education, 
Culture, Sports, Science and Technology. 
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Fig. 1. — Photoabsorption cross-section of Carbon, Oxygen, and Iron. The open circles 
in dicate those obtained fror n the NIST database. The sohd and dash ed hnes show those 
by IVerner et al. I (119961 ) and iBalucinska-Church &: McCammon I (119921 ). respectively. Their 
cross-sections are referred to as vern and heme in XSPEC. 




Fig. 2. — A cross section view of the torus structure in our simulation. The primary X-ray 
source is located at the center of an obscuring torus. The torus structure is characterized by 
the half-opening angle ^oa, the inclination angle of the torus from an observer 9i, the column 
density A^h along the equatorial plane, and the ratio of Tin to rout- A simulated spectrum 
is separated into three components: one direct component and two reflection components. 
The two reflection components are referred as reflection component 1 and 2. The reflection 
component 2 consists of reflection light from the visible inner wall of the torus. The reflection 
component 1 consists of the rest of the reflection light. 
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Energy (eV) 

Fig. 3. — An example of our simulated AGN spectrum with A^h = 10^^ cm"^, = 40°, and 
= 45°. The simulated spectrum is divided into three components : the direct component, 
the reflection component 1, and the reflection component 2 as shown in this figure. The 
dashed line displays the intrinsic source spectrum. 



g Direct component g Reflection component 1 g Reflection component 2 




Energy (eV) Energy (eV) Energy (eV) 



Fig. 4. — A^H dependence of the three components for = 40° and 6i = 45°. The column 
density has values of A^h = 5x10^3, IQ^\ 2x10"^^, 5x10^4 and lO^^ cm'^. 
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Fig. 5. — Comparison between the simulated direct component and the cut-off power law 
model affected by both photoelectric absorption and Compton scattering. 



o N„=10^ (cm-2), inclination angie = 89-91 (deg) o N„=102'' (cm-^), inciination angie = haif-opening angie + 1 (deg) 




1000 10^ 10= 1000 10' 105 



Energy (eV) Energy (eV) 

Fig. 6. — ^oa-dependence of the reflection components 1 and 2. For the simulations of the 
reflection component 1, A^h and 9i were fixed at 10^"^ cm~^, and 90°, respectively. For the 
simulations of the reflection component 2, 9i was flxed on ^oa+l°- 
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Fig. 7. — Total counts of the reflection components 1 and 2 as a function of ^oa- Open 
circles and open squares represent the reflection component 1, and the reflection component 
2, respectively. 
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Fig. 8. — The Qoa dependence of the reflection component 1 (left) and 2 (right) for various 
column densities, A^h=10^^, 10^^, and 10^^ cm~^. The left panel shows the normahzed count 
rate divided by that simulated in the spherical distributed matter. The solid lines in the left 
and right panels show the covering factor of the torus, cos^^oaj and the combination of the 
covering factor and the apparent size of the visible inner wall, cos6'oa cos(26'oa ~ 7'"/2). 
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o N„=102'' (cm-2), half-opening angle= 10 (deg) o N„=10^' (cm-^), half-opening angle= 10 (deg) 




1000 10' 105 1000 10" 10= 



Energy (eV) Energy (eV) 

Fig. 9. — The ^i- dependence of the reflection components 1 (left) and 2 (right). In these 
simulations, A^h and 9oa. were set at 10^^ cm~^ and 10°, respectively. 
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Fig. 10. — Total counts of the reflection components 1 and 2 in the 1 keV-300 keV band as 
a function of 6i. Open circles and squares show the total counts of the reflection component 
1 and 2 at ^^oa=10°, respectively. We also show the total counts for ^oa= 30° by closed circles 
and squares. A curve of cos 9i is plotted in the dotted line, in order to show a 6'i-dependence 
of the projected area. 
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Fig. 12. — The fraction of the absorbed luminosity with respect to the intrinsic source 
luminosity. The fraction was estimated for Nh=10^'\ lO^'l 10^^ and 10^^ cm . We also 
plot a curve of cos 6'oa in the dashed line, in order to display the covering factor of the dusty 
torus. 
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Fig. 13. — Equivalent width of the iron Ka line to the total continuum emission, summed 
over the direct and reflection components. The left panel shows the A^h dependence of the 
equivalent width for ^^oa=30°, and the right panel shows the dependence on 6oa for iVH=10^^, 
10^^, 10^^, and 10^^ cm^^ from the bottom to the top (or from thin to thick lines). The 
open circles, squares, and triangles indicate the equivalent width in the ranges ^^i=0-18°, 
60-63°, and 81-84°, respectively. We plot the equivalent width for the spherical distribution 
at 6'oa=0, instead of using our simple torus model with ^oa=0. 




Fig. 14. — Equivalent width of the iron Ka line to the reflection components 1+2. Symbols 
are the same as in Figure 12. 
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Fig. 15. — The normalization of the direct and reflection components as a function of ^oa- 
^i=6'oa+l° is assumed. The closed and open circles display their normalizations as shown in 
the figure. The solid hne indicates a curve proportional to l/cos6'oa- 




Fig. 16. — The wide-band spectrum of Mrk 3 observed with Suzaku fitted with our model 
(left), and the unfolded X-ray spectrum (right). The spectra in the left panel are obtained 
by the XIS-FI (black), XIS-BI (red), and HXD-PIN (green). The energy range around the 
Si-K edge (1.875-1.845 keV) is ignored for spectral fitting. The right panel shows continuum 
emission from PLl, reflection 1, reflection 2 and heavily absorbed PL2. 
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Table 1: Parameters and grids of tabic iriodcl 



Parameters 


Parameter grids 


Photon index 


1.5, 1.9, 2.5 


Nn (xl022 cm-2) 


1, 5, 10, 50, 100, 200, 300, 500, 700, 1000 


half-opening angle (°) 


0, 10, 20, 30, 40, 50, 60, 70 


inclination angle (°)" 


1- 89 in steps of 2 



"We selected X-rays within the inclination angle± 1°. 



Table 2: x' •-^li tl^f- parameter grid of O^j.^ and 6i. 



di — ^oa 










^oa 












30° 


35° 


40° 


45° 


50° 


55° 


60° 


65° 


70° 


-2 


704.6 


706.7 


706.2 


702.1 


701.9 


700.8 


706.6 


703.3 


706.3 


-1 


701.6 


705.9 


705.7 


702.8 


700.2 


703.2 


708.9 


704.0 


706.2 





701.7 


705.3 


706.2 


704.1 


702.9 


707.3 


707.0 


706.7 


707.6 


+ 1 


705.1 


710.2 


711.5 


708.8 


711.8 


710.2 


713.4 


713.6 


723.6 


+2 


724.9 


736.4 


724.4 


734.1 


729.3 


737.6 


737.7 


752.0 


753.6 



Table 3: Best-fit model parameters. 



Photon Index A^hi 

(xl024 cm-2) 


T direct 
-'^2-10 

( erg s-i ) 


-^Vh2 ^2*^10 ^oa 

(xl024cm-2) (ergs-i) (°) 


9i 

n 


xV(d.o.f.) 


1.82 1.1 


2.1 xlO'*^ 


1.5 5.1x10^3 50(fixed) 


51 (fixed) 


726/613 



Note. — A^Hi was linked to A'h2. 
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